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We use a version of the density functional theory to study the solvation force between two plates modified with 
a tethered layer of chains. The chains are built of tangentially jointed charged spherical segments. The plates are 
immersed in an electrolyte solution that involves cations, anions and solvent molecules. The latter molecules 
are modelled as hard spheres. We study the dependence of the solvation force and the structure of chains and 
of solute molecules on the grafting density, length of chains, architecture of the chains and on concentration of 
the solute. 
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1. Introduction 

Layers of charged polymers tethered onto solid surfaces have attracted great attention for the recent 
years due to their biological importance and numerous practical applications OjJ. The behavior of charged 
olymers at interfaces is not only more complex as compared to charged polymers in bulk solution la- 
" but is also more complex than the behavior of systems involving uncharged tethered polymers |6|-|8j]. 
Generally, these systems contain polymers carrying positive and/or negative charges, surface charges, 
mobile ions in solution and solvent molecules. All these components make the properties of a system 
intricately dependent on electrostatic interactions jil-tlill. 

A special type of charged polymers are polyampholytes, i.e., polymers consisting of positively and 
negatively charged segments and, in some cases, also neutral segments. Differently charged, as well as 
neutral segments, can be distributed randomly or in a certain sequence ll4T-[l7ll. A number of synthetic 
and natural (e.g., proteins) polymers can be classified as polyampholytes ll8U20ll A special example of 
such systems are the so-called polyzwitterionic polymers I21I1 . 

Due to their unique features, polyampholytes have a wealth of important and practical applications, 
e.g., as dispersing additives, gelling agents, rheology modifiers, etc. I18II . Nevertheless, polyampholytes 
have received much less attention, particularly in theoretical studies, as compared to neutral polymers. 
Consequently, our understanding of polyampholytes, especially tethered polyampholytes, is far from be- 
ing complete. 

In the case of uncharged solid surfaces and tethered chains whose segments bear the same charges, 
the electrostatic repulsion forces can be strong even at low grafting densities. Consequently, under such 
conditions, the chains are stretched and the system enters the brush regime [17]. In the case of tethered 
polyampholytes, the pinned layer behaves like a polyelectrolyte when the chains are globally charged 
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and the resultant single chain charge is significant. However, for zero or for low net charge, the polyam- 
pholyte can coil and the tethered layer collapses. Low density tethered polyampholytes exhibit the highest 
deformation amplitude among the stretched and the collapsed states |22]|23]]. The structure and thermo- 
dynamic properties of tethere d polye lectrolytes have been investigated using both theoretical 
and computer simulation Ilall7ll30l-l35ll methods. Density functional methods are among the most suc- 
cessful theoretical approaches I36U41I1 . 

Interaction forces between surfaces and colloidal stability are closely related. A stable colloidal sus- 
pension is characterized by repulsive forces between the colloidal particles, while attractive forces lead 
to an unstable suspension. One might expect that pinned polyelectrolyte layers induce repulsive forces 
and thus it leads to suspension stabilization. However, it has been also observed that polyelectrolytes 
may destabilize colloidal suspensions [42]. The problem of evaluation of the solvation forces between the 
surfaces modified with tethered charged chains and their dependence on the properties of a solvent and 
on the properties of the chains themselves was discussed in numerous papers 1341 143114911 . 

The purpose of this work is to use the density functional theory to calculate the force between two 
plates covered with tethered charged brushes. The approach used by us is an extension of the theory 
outlined in reference fiill . We will consider both the chains with the segments having identical charges 
(i.e., typical polyelectrolytes) and polyamphylic molecules with no net charge per chain. In the latter case, 
we wonder how the distribution of the charges along the chain (i.e., the chain architecture) effects the 
solvation force and the structure of a system. 



2. Theory 

We consider two identical surfaces, lower and upper, at a distance of H apart. Each surface is covered 
with charged tethered chains (P). The chains are built of M tangentially jointed charged hard spheres of 
the same diameter, a. The charge of the segment j, j = 1,...,M is zj P) e, where e is the magnitude of an 
elementary charge. The chain connectivity is ensured by imposing a bonding potential I36H39I1 



M 



Vb(R)= £wb(|r J+ i-r y |), 

1=1 



(1) 



where R = (n,.. . ,rjvf) is the vector describing the positions of all the segments. For tangentially bonded 
chains, the potential Vj, satisfies the relation 



exp 



-Vb(R) 
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The first segment of each chain, j = 1, is pinned at one of the surfaces. The surface bonding potential 
has the form 



exp 



kT 



■■ ^S(z'), 



(3) 



where ^ is a constant and z' - z-a/2 for lower and z' - z - (H-a/2) for the upper surface. The inter- 
action of all the remaining segments of a chain with the walls is described by the potential 



vj p \z)=vf\z) + v^z), 



(4) 
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where v . (z), j = 2, . . . ,M, is the hard-wall potential 



oo, for z<a/2 or z> H - a 12, 
0, otherwise, 
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and v^(z) = V j>el (z) + V}, e i(H - z), where 
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is the Coulomb potential. In the above, Qe is the surface charge density of each wall and e is the dielectric 
constant. According to equation {6j, the surface charge density on each wall is assumed to be the same. 
The total chain-surface potential is then V {P \z x , . ..,Zj) = Y.f =l vf\zj) + T.f =l vf^{zj). 

The modified surfaces are in contact with a three-component fluid. The indices a = S, C and A abbre- 
viate solvent, cation and anion species, respectively. The particles of the second and of the third species 
(C, A) are charged hard spheres as well, for the sake of simplicity, of a diameter a, and carry the charges 
eZy and eZj\ In other words, we consider the so-called solvent primitive model of an electrolyte. This is 
the simplest model of electrolyte solutions that takes into account the presence of non-zero size solvent 
molecules. The polar nature of a solvent is taken into account by retaining the dielectric constant in the 
Coulomb interaction between the ions. Note that the solvent primitive model was used by Henderson 
et al. to describe the forces between macroscopic particles in an electrolyte and to model electrolyte- 
membrane systems ji^lijl. Other examples of applications of this model can be found in references IS2| — 

m 

Despite its extreme simplicity, the solvent primitive model has at least two virtues. It may be as rea- 
sonable a model of a solvent as can be treated conveniently by density functional theory and it does 
recognize that solvent molecules occupy space. 

In this work we restrict our attention to the case of monovalent ions, i.e., Zj = |Zp | = 1 and \Zj\ = 1 
for j - 1,...,M. The solvent S molecules, however, are uncharged hard spheres of the same diameter, 
a. Thus, = 0. The interactions of the fluid species with the surface are given by equations {D~GD, m 
which zj p) should be replaced by Z® . 

The interactions between all the spherical species (i.e., between the molecules S, A and C and the 
consecutive segments of the chains P) are given by 

, oo, r <a, 

uT{r) = \ e^z^z" (7) 



r> a, 



where 77, a = S, C,A,P. 

The confined system is in equilibrium with a three component bulk mixture of the components C, A 
and S. The bulk system is at the same temperature and at the same chemical potentials ji^, rj — A,C,S as 
the confined system. The bulk densities of particular species are pjj. 

In order to proceed, let us introduce the notation, p (P) (R) and p tr,) (z), 77 = S,C,A for the density dis- 
tribution of chains and of fluid species, respectively. The theory is constructed in terms of the density of 
particular segments of chains, p (P) (z), and the total segment density of chains, pf\z). These densities 



are introduced via commonly used relations I36|439[|41H 

M M f 

pT w = E p l s P i w = E dR5 ( r ; - r )p (P) (»)■ (8) 

In the system under study, all the local densities are only functions of the distance z from the lower 
surface that is set at z = 0. 

The system is studied in a grand canonical ensemble with the constraint of constancy of the number 
of tethered chain molecules, i.e., 

p c = J dzp^Hz), (9) 

where p c is the total number of tethered chain molecules per area of the surface. Since the confined 
system comprises two surfaces, the surface density of the chains is Rq = p c /2. 

The equilibrium density profiles are obtained by minimizing the thermodynamic potential 

dRp (P) (R) E v [P) {Zj)+ E Av[v {r >\z)p^\z)-p 7] \+\ Avq{z)W{z). (10) 

j = l 77=S,A,C>' J 

In the above, F is free energy functional and q{z) is charge density 

— = t^p P fhz) + E zfp^Hz). (id 

e j = l t;=S,A,C 
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The electrostatic Wiz) satisfies the Poisson equation 



9 An 
V l *¥{z) = q(z). 



(12) 



A solution of the differential equation (12) for a fluid confined between two walls is given in refer- 
ence (il|]. It requires denomination of the value of electrostatic potential at the wall, = ¥(0) = ¥(//). 
From the electro-neutrality condition of the system, it follows that 



/ 



Q+ / &zq{z) = 0. 



(13) 



The principal task in applying the density functional theory is to derive an expression for the Helmholtz 
energy as a functional of the density profiles. Following previous works , the Helmholtz en- 

ergy is divided into an ideal term that depends on the bond potentials and the architecture of polymers, 
and the excesses arising from interactions of various forms. The latter terms are responsible for the 
thermodynamic nonideality. Within the framework of our model, the Helmholtz energy functional, F, is 
decomposed into the following contributions, 



F = F [d + F ex , 
where Fid is the ideal contribution, 



Fex- F P + F hs + F ei , 



^ = J-JdRp (p >(R)V b (R) + jdRp^(R){ln[p< p) (R)]-l}. 



(14) 



(15) 



The volume exclusion (the hard-sphere) term, Fhs. is calculated according to the Fundamental Measure 
Theory, cf. references (^.[i^liill. 
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where «/, i = 0, 1,2,3 and nyj, j — 1,2 are, respectively, scalar and vector total weighted densities. The 
total weighted densities are sums of the weighted densities of individual species. For example, 



n^nF + nP + nP + nP, 



(17) 



and Since the relevant equations defining the weighted densities have been already 

presented in numerous works ItiHil liiU^ . we have omitted them here. 

The contribution _F e i, arising from the coupling between electrostatic and hard-sphere interactions is 
written down employing the approach described in detail in references (i^li^l We have 



Fa = a>ei(z)dz, 



where 



®elU) = - — 



M 
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r r 3 
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1 + oT 3n 2 ' 

In the above T = ( V1 + 2k<j- l) 12a, T* is the "electrostatic" reduced temperature, T* = kTea/e 2 , 



K 2 = 
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(18) 



(19) 



(20) 



The quantities pS) and pW a re the "ref erence electrostatic system averaged densities", which are calcu- 



lated according to references I39LI57L I5 
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Finally, the excess free energy functional due to the intra-chain correlation is given by 



FplkT ■■ 



1 



where ^ = l-n^-V2^ 



1-M 

>2 



M-l n M c 
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(21) 



/ ( W 2 P) ) anc * yj * s tne contact value of the cavity correlation function. 
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(22) 



We are aware of the fact that in our treatment, the calculation of the values of T is performed without 
discriminating the ions belonging to chains and free ions in the confined solution. In the light of the 
works such ions should be distinguished. Therefore, the above expressions fL9l-j22), 

and equation (22) in particular, are approximations. Note that in the case of bulk fluids, the sums in 
equation (21) are identical. 

At equilibrium the density profiles minimize the thermodynamic potential &, i.e., 



81V 



<5p( p) (R) ~ 8pW(r) 
This condition leads to the equations 



= 0, ?7 = S,A,C. 
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(23) 
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J 3 j 



rCP)\ 



SFe, 



+ v w (z) + eZ w '¥{z), 77 = S,A,C, 



(26) 



and where the constant C is calculated from the normalization condition (9). The multidimensional den- 
sity profile equation J24t can be then reduced to the equations for local densities of consecutive segments, 
p[ p ) (z/) using the method described in fiill . 



The solvation force per unit area is calculated from 



f s _ d<3/(H)/AkT p 
kf ~ dH kT' 



(27) 



where A is the surface area and p is the pressure of the bulk deference fluid involving the components S, 
A and C. 



3. Results and discussion 

The systems in question are characterized by numerous parameters. In order to reduce their number, 
all the calculations have been performed assuming constant total bulk density of the fluid, pj^ = p^a 3 = 
0.7, where pb = p b A) + p b C) + p b S) . This value is close to the density of water solutions at standard temper- 
ature and pressure (STP) conditions. The composition of the bulk fluid, i.e., the bulk mole fraction of the 
1:1 electrolyte, x = (p^' + p b C ')/pb. was varied. Also, the value of the reduced electrostatic temperature 
was kept constant and equal to T* = 0.2. This temperature is lower than the temperature usually used in 
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computer simulations of nonuniform electrolytes, but our calculations have been intentionally carried 
out at a lower temperature in order to increase the role of electrostatic interactions. The parameters that 
characterize the brush are the number of segments and the surface brush density, R^, = Rqcf 3 - Numerous 
calculations were carried out assuming that the total charge of the brush is zero, i.e., that the number of 
positively and negatively charged segments is the same, M+ = M_ = MI2. However, the distributions of 
"+" and "-" charges along the chain were different. We have also studied the cases of all the segments of 
the brush bearing the same charges. Moreover, in some cases, we carried out calculations for uncharged 
chains and for the chains with non-zero resultant charge, lower than Me. In all cases, the electrostatic 
potential at the wall, To. was set to be zero. 

We introduce the following codes to distinguish the systems under study. The symbol m\ + n\ - rri2 + 
ri2-,... abbreviates a chain whose first mi segments are positively charged, the next n\ - negatively 
charged, etc. When the sequence of the charges along the chain is repeated, we use parentheses to group 
the repeating units, for example the symbol 5(1 + 1-) means that the chain is built of 10 segments alter- 
nately charged with + and -, and whose first (pinned) segment is positively charged. 




H 



Figure 1. (Color online) Part (a). The solvation force vs. H* = HI a. All the results, except for those marked 
by a dotted line that were obtained for uncharged chains (M = 10), are for 10+ chains. Solid line is for 
R* = 0.0025, dash-dotted line - for R* = 0.1 and dashed line - for R* = 0.15. Part (b). Examples of the 
local densities of the chains 10+ and of the fluid. Upper panels are for R^, = 0.0025 and H* = 3 and 2, 
while lower panels - for R^, = 0.1 and H* = 10. Left upper panel: solvent density profiles and the profiles 
of charges, due to ionic species. Right upper panel: the total segment density profiles. The calculations 
are for H* = 3 (dashed lines and open symbols) and for H* = 2 (solid lines and filled symbols). Lower 
left panel: the total segment density profiles of charged (solid line) and of uncharged chains (dashed line). 
Right lower panel: the profiles of positive (dashed line), negative (dash-dotted line) ions and solvent (solid 
line, the latter profile is multiplied by 10) for 10+ brush. For uncharged brush, the ionic species of profiles 
are identical (dash-doubly dotted line). The solvent profile is given by a dotted line. The latter profiles are 
multiplied by 10. The bulk mole fraction of the electrolyte is x = 0.1. 

In figure[T](a), we show how the solvation force depends on the surface density of the grafted chains 
10+, whose all segments bear the unit positive charges. For a comparison, we also display here the result 
for uncharged chains built of M- 10 segments. For a low surface density of the chains, = 0.0025, the 
course of the solvation force vs. H* reminds us the course for the pore with non-modified walls The 
solvation force diverges at low wall-to-wall separations H* - HI a and exhibits oscillations correspond- 
ing to attractive and repulsive forces that are well-pronounced at small values of H* and decay to zero 
at larger plate-to-plate separations. The maxima of the solvation force correspond to the development of 
consecutive layers of the solvent density profiles, see figure[T](b). The role of chains and of electrostatic 
interactions in particular is small and the solvation force is almost entirely determined by the solvent 
packing effects. 

When the surface density of chains increases, the dependence of the solvation force on H* quanti- 
tatively changes. For charged chains (10+), the solvation force becomes repulsive for all values of H* . 
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Two factors play a significant role here. The first one is the volume exclusion effect: with an increasing 
Rc, the segments of the chains occupy more and more space, and further compression becomes difficult 
for small H* . The second factor is electrostatic repulsion between the segments of chains pinned to the 
opposite walls. For R~ = 0.1, the second effect becomes very important. Indeed, the solvation force for 
charged and uncharged brush of the surface density 0.1 is very different. For uncharged chains [dotted 
line in figure[T](a)], the solvation force becomes even slightly attractive (e.g., for H* « 4). For H* s= 3 the 
solvation force becomes repulsive and grows rapidly. Electrostatic forces between the segments cause 
the appearance of strong repulsion between two walls at larger distances H*. The above results are in 
quantitative agreement with experiments fiilliill . 

Figure [1] (b) shows examples of a structure within some selected pores. Upper panels are for i?£ = 
0.0025 and for the chains 10+. Left upper panel displays the solvent density profiles (solid lines) and the 
negative total charge profiles (symbols) for two pores with H* = 3 and H* = 4. For H* = 4, three well- 
developed layers of the solvent appear within the pore. This structure is resistant to compression and a 
local maximum appears on the plot of the solvation force. When H* decreases down to 3, the inner layer 
is "squeezed out" the pore, and the next structure develops which is resistant to further compression. The 
total segment density profiles [the right upper panel of figure[l](b)] behave similarly to the profiles of the 
solvent. 

The average density of a fluid within the pore is much lower than in the bulk reference system (we 
recall that the latter equals 0.7). A part of the space inside the pore is occupied by the segments of chains 
[see the upper right panel of figure Q](b)], but electrostatic interactions cause a lower average density in 
the confined system that is in chemical equilibrium with the bulk system. The upper left panel shows 
the negative charge density profile of ionic species of the fluid, -q'*{z) - -q'(z)a 3 le, where q'{z) = 



Ia=A,Cp W U)- 

Despite the fact that the electrostatic potential at the pore walls is zero, = 0, an effective nega- 
tive charge on the pore walls is generated. This is the result of the charges on the segments. Unlike the 
pores with non-modified walls (i?£ = 0), the zero surface charge does not correspond to the zero of the 
electrostatic potential at the wall. It would be of interest to determine the dependence of the so-called 
"potential of zero charge", PZC, on the pinned chain characteristics, but this problem is out of the scope 
of the current research. 

Lower panels in figure[T](b) show the structure inside the pore of H* - 10 wide. The grafting density 
if high, i?£ = 0.1. Due to the symmetry, only one half of the profiles is shown. We compare here the 
situation for uncharged chains built of M = 10 segments and for the chains 10+. Uncharged tethered 
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Figure 2. (Color online) The solvation force vs. H* = HI a. Solid line is for the 10+ chains, dashed line - 
for the chains 5(1 + dotted line -for 4(2 + 2-) 1 + 1- and dash-dotted line -for 2(5 + 5-). Part a is for 
Rt = 0.01 and part b - for Rt - 0.1. In all cases the bulk mole fraction of the electrolyte is x = 0.1. 
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chains are much more coiled. Indeed, their density is low at the pore center. The stretching of charged 
chains is caused by electrostatic repulsion between the segments. Positively charged segments attract 
negative ions, which, in turn, cause the accumulation of positive ions. However, the concentration of 
solvent molecules inside a pore is very low, because there is simply "no room" for them. In the case 
of uncharged chains, the concentration of all components of the fluid is very low inside the pore. This 
is because there is no attraction between the segments and ions, and the volume exclusion effects are 
responsible for a very low confinement of all the fluid species. 

Here we consider the cases of the total charge of a single chain equal to 0. We studied the following 
distributions of the charges along the chain: (a) 5(1+ (b) 4(2 + 2-) 1 + 1- and (c) 2(5 + 5-). The brush 
10+ is treated as a reference. Figure[2](a) shows the solvation forces for four systems defined above. The 
surface grafting density was - 0.01. Surprisingly, despite quite big differences between the considered 
chains (chains with the total net charge, but with different distribution of charges vs. the chain with 
the total charge equal to lOe), the differences between the solvation force for particular systems are 
rather small. The solvation force for the chain 10+ is more repulsive. The maxima of / s are higher and 
the minima are shallower than for all the remaining cases, as expected from the consideration of the 
role of electrostatic forces. Quite unexpected, however, is a very small effect of different distributions of 
the charges along the chain on the solvation force. Indeed, the results for the systems (a)-(c) are hardly 
distinguishable on the figure scale. At a higher grafting density, = 0.1, the solvation forces for the 
systems (a)-(c) are still very similar. However, the solvation force for the chains 10+ is very different 
from all the remaining curves and shows that an electrostatic repulsion occurs at larger wall-to-wall 
separations, cf. figure[2](b). The pinned chains (a)-(c) can assume coiled configurations at both pore walls, 
whereas for the chain 10+, electrostatic repulsion between segments inhibits its coiling. 




z z z 



Figure 3. (Color online) Part (a). Left panel shows the total segment density profiles for the systems in 
figure[2](a). Abbreviations are the same as in figure[2](a). Right hand panel shows the profiles of ions for 
two systems shown in figure[2](a), namely for 10+ (solid lines) and 5(1 + 1-) (dashed lines). Lines without 
symbols are for anions, a = A, lines with symbols are for cations, a = C. Part (b). Left hand panel shows 
the solvent density profiles for grafted chains 10+ (solid line) and 5(1 + 1-) (dashed lines), while the right 
hand panel shows the total charge density profiles, q*{z) = q{z)a s e. Abbreviations are the same as in 
figure[2](a). Due to the symmetry, only one half of the profiles is displayed. The pore width is H* = 10, 
and the bulk mole fraction of an electrolyte is x = 0.1. 

Examples of a structure for the systems shown in figure [2] (a) are shown in figure [3] Although the 
width of a pore is quite large, H* = 10, there is still a significant overlap of the chains pinned at the 
opposite walls, especially for the 10+ chains. Despite small differences between the solvation forces [fig- 
ure [2] (a)], the differences between the total segment density profiles [right panel in figure [3] (a)] for the 
cases (a)-(c) and 10+ are pronounced. Although the latter chains are more stretched than all the remain- 
ing chains, the total segment density profile for the chains 10+ is rather low at the pore center. Due to a 
low surface density, there is enough space for the segments to assume a slightly tilted or parallel to the 
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wall configuration. Right panel in figure[3](a) shows the profiles of ions, a — A,C. A significant difference 
between the profiles of anions and cations is observed for the chains 10+ due to electrostatic segment- 
anion attraction. For electro-neutral chains, the differences between the profiles of anions and cations 
are small. However, the profiles of ionic species are not identical. This indicates that there should be a 
small extra charge on the pore walls. In other words, the PZC for electro-neutral chains is not at 'Prj = 0. 
An interesting question is how PZC depends on the architecture of the chains, but this problem is out of 
scope of the current research. 

Figure [3] (b) shows examples of the solvent density profiles (right hand panel) and the charge den- 
sity profiles (left hand panel). Distribution of the solvent particles inside the pore depends very little on 
the charges of the segments. The solvent particles "do not feel" the charges, and the volume exclusion 
effects almost completely determine the structure of the confined solvent. The total charge distribution 
is sensitive to the architecture of the grafted chains. 

Similar investigations of the structure have also been carried out for a higher grafting density, Rl - 
0.1. The solvation forces for those systems are shown in figure[3](b). Figure|4](a) compares the total seg- 
ment density profiles for the systems 10+ and for the systems (a)-(c) (left hand panel), as well as the 
density profiles of fluid components for the systems 10+ and 5(1 + 1-) (right hand panel). The pore is 
narrow, H* - 5, so the segments of the chains are compressed and form layered structures. The differ- 
ences between the total segment densities for different systems are not big. Because of a small pore width, 
there is not much room for different rearrangements of the segments and the observed layered structure 
permits to effectively pack the segments inside the pore. 





Figure 4. (Color online) Part (a). Left panel shows the total segment density profiles for the systems in 
figure[2](b). Abbreviations are the same as in figure[2](b). Right panel shows the profiles of ions for two 
systems from figure[2](b), namely for 10+ (solid lines) and 5(1 + 1-) (dashed lines). Lines without symbols 
are for anions a = A, with symbols - for cations, a = C. The latter profile was multiplied by 100 to make it 
visible. We also shown here the solvent profile for the system 5(1 + 1-) (dotted line, a = S). Part (b). The 
profiles p^' (z) (left panel) and the total charge profiles (right panel) for the systems from figure[2](b). The 
meaning of the lines is the same as in figure[2](b). The pore width is H* = 5 and the bulk mole fraction of 
the electrolyte is x = 0.1. 



We have also inspected the density profiles of individual segments, but for the sake of brevity only 
the profiles of the middle segments, i - 5, are shown in the left hand panel of figure|4](b). The differences 
between these profiles are significant. The performed analysis of the shape of the functions p^'(z), i - 
1,2, . . . , 10 leads to the conclusion that in the case of 10+, the layers of chains pinned at the opposite pore 
walls penetrate into one another. However, in the remaining (a)-(c) cases, the chains "turn back" at the 
middle of the pore and the inter-penetration of the layers attached to the opposite walls is much weaker. 

In the case of the pore H* = 5 modified with 10+ chains, the confined fluid contains almost solely 
anions (figure [4] (a), left hand panel). The concentration of cations is extremely low; in order to make the 
cations profile visible on the figure scale, we multiplied it by 100. Similarly the density of the solvent 
inside the pore is also extremely low. In the case of a pore modified with 5(1 + 1-) brush, the anions 
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accumulate at the pore walls, while in the pore interior cations prevail. Moreover, for 5(1 + 1-) chains, 
we also observe a concentration of solvent molecules inside the pore. 

The shape of total charge distributions inside the pore exhibits quite large local changes. In all the 
cases studied, there appears a very high peak at z* - 0.5 on the curves q{z) (we have cut the height 
of this peak on left hand panel of figure [4] (b) in order to make the figure readable). The height of this 
peak almost exactly corresponds to the amount of pinned charged segments (in all the cases, the pinned 
segments are positively charged). For a system 10+, the function q{z) is positive within the region of 
0.5 < z < 1.5, while for the same region it is negative for all the remaining systems in question. Similarly 
to the previously considered cases, the integrals of q{z) over the entire pore are different from zero, i.e., 
there should be extra charges at the pore walls in order to ensure the total electro-neutrality of the entire 
system. In other words, the PZC does not correspond to the zero value of electrostatic potential at the 
wall. 




Figure 5. (Color online) Solvation force for the brushes 8+ (solid line), 4(1 + 1-) -dashed line, 2(2 + 2-) - 
dots and 4 + 4- and dash-dotted line. Part (a) is for i?* = 0.05, part (b) is for i?* = 0.025. The bulk mole 
fraction of the electrolyte is x = 0.1. 

Now we consider the chains built of M — 8 segments. Similarly, we study the case when all the seg- 
ments are positively charged, 8+, as well as three electro-neutral chains: (a) 4(1 + 1-); (b) 2(2 + 2-) and 
(c) 4 + 4- . The calculations of the solvation force have been carried out for several grafting densities R^, 
but only some selected results are shown in figure [5] In general, all the features observed for the chains 
of 10 segments also appear here. When the grafting density is not too high, the solvation force oscillates 
around zero. For higher values of the solvation force becomes entirely repulsive, especially in the 
case of chains 8+. For electro-neutral chains, the solvation force only weakly depends on distribution of 
the charges along the chain. In the latter case, the difference between the solvation force for different 
architectures of the brush slightly increases with an increase of the grafting density [cf. figure [5] (a) and 
figure[5](b)]. Since the changes in the structure of a confined fluid with distribution of charges along the 
chain are quantitatively similar to M — 10 cases, the relevant figures have been omitted. 

Figure [6] shows the dependence of the solvation force on the bulk mole fraction of an electrolyte. The 
calculations are for two types of grafted chains, namely for 4(1 + 1-) [part (a)] and for 8+ [part (b)]. In both 
cases, the effect of the bulk concentration of ions on the solvation force is small. It is slightly more pro- 
nounced for 8+ chains, especially for 1 1.5 < H* < 2. In other words, the most important factors effecting 
the solvation force are the characteristics of the chain, i.e., its length, grafting density and architecture. 

The results presented above have shown that one of the important factors effecting the behavior 
of the solvation force is the resultant charge of a brush. Therefore, in figure [7] we present the results 
for the chains of 8 segments with the resultant charge equal to 2e, 4e, 6e and 8e, respectively. Now, the 
differences between the solvation forces for a particular system are big. For a larger resultant charge a 
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H H 

Figure 6. (Color online) The dependence of the solvation force for the systems 4(1 + 1-) [part (a)] and 8+ 
[part (b)] on the bulk mole fraction of the electrolyte. Solid line is for x = 0.2, while dashed line is for 
x = 0.05. The grafting density is R* = 0.025 in part (a) and 0.5 in part (b). 




Figure 7. (Color online) Solvation force for the systems 5 + 3- (solid line), 6 + 2- (dashed line), 7+1- 
(dotted line) and 8+ (dash-dotted line). The grafting density is Rt - 0.05 and the bulk electrolyte mole 
fraction is x = 0.1. 



strong repulsion occurs at larger wall-to-wall separations. This is obviously due to electrostatic repulsions 
between the chains tethered at the opposite walls. When the resultant charge of a chain is 2e, we still 
observe effective attractions between two pore walls at H* a 2.5 and for H* a 3.5. For the resultant 
charge of 4e, there appears a marginally attractive force between two walls at z « 2.5. However, for larger 
resultant charges, the solvation force is repulsive for all H* , although it exhibits several local minima and 
maxima. When the surface grafting density increases, the local extrema become weaker, and the curve 
describing the dependence of the solvation force on H* becomes smoother [cf. also figure [2] (a)]. 
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4. Summary 

In this work we applied the density functional theory to the study of the solvation forces between the 
walls covered with tethered layers of charged chains. We considered the chains possessing different net 
(resultant) charges and investigated the cases when all the segments had the same charges, as well as 
polyampholytes with zero and non-zero resultant charge per chain. The results of calculations showed 
that one of the factors that effects the solvation force very much is the resultant charge of the chain. 
However, the distribution of the charges along the chain effects the solvation force rather little. Another 
factor that quantitatively changes the solvation force is the grafting density. On the other hand, the con- 
centration of the solution plays rather a minor role. Of course, the above results are valid for a specific 
model studied in this work. Moreover, we considered rather short chains and thus cannot exclude that 
our conclusions regarding the tethered layers built of chains involving hundreds of segments would be 
different. 
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Ci/ma cosibBaTaLjiT Mi>i< ma pa mm npi/ieflHam/ix no/iie/iei<Tpo/iiTiB. 
niflxifl 3a flonoMoroK) <(>yHKi^ioHa/iy rycTi/iHi/i 

0. ni3icR A. naTpi/iKee^, C. Coko/iobckR fl. I/ibH 1/114 bKi/iiP 

iHCTHTyT xi m i Y yrAM, 04360 MexiKO, MeKcnKa 

Biflfli/] MOfle/iiOBaHHfl <()i3WKO-xiMiHHHx npoijeciB, yHiBepcuTeT MapiT Kropi-CioioflOBCbKoT, 
20031 J]k>6ji\h, no/ibina 

iHCTHTyT cf>i3HKH KOHfleHCOBaHnx CHdeM HAH YKpaiHn, syn. CB£HL4ii4bKoro, 1, 7901 1, /lbBiB, YKpai'Ha 

3a flonoMororo BapiaHTy Teopn (fiyHKLjiOHa/iy rycTWHH bmbhcho cn/iy co/ibBaTaijM Mix flBOMa n/iacrnHaMn, no- 
BepxHi akhx MOflW(()iKOBaHO LuapoM npweflHaHwx no/iiMepmix zianmtoxKiB. OdaHHi ccfiopMOBaHO y Bwnnafli 3'efl- 
HaHnx flOTH4Hwx oflHH flo oflHoro 3apflflxeHwx c(()epw4Hnx cerMeHTiB. n/iac™™ 3aHypeHo b po3HWH e/ieiapo- 
n'ny, rkww MicTHTb KaTiomi, aHioHn Ta MO/ieKy/in po3HWHHHKa (si Ki MOfle/iKDKDTbCfl hk TBepfli ccfiepn). Mm koh- 
L4eHTpyeMOCb Ha 3a/iexHOC"ri cn/in co/ibBaTaijM Ta CTpyicrypn hk no/iiMepi-inx /iaHi4K)XKiB TaK i po34HHHHKa Bifl 
rycTWHW npneflHaHHfl /iam;K)XKiB, i'x aobxhhh, apxiTeiaypn Ta KOHL^eHTpai^iT po34HHeHO'i penoBUHn. 

IOi K) h o b i c/iOBa: np\Aep,naH\ eneKiporinw, cmna conbBajau,i'i, ap,cop6u,\n, reopin cpyHKLfiOHayiy rycrnHH 
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